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Chiral shifts in heavy-light mesons 

A.M.Badalian, Yu.A.Simonov, M.A.Trusov a 
a ITEP, Moscow, Russia 

The mass shifts of the P-wave D s and B a mesons due to coupling to DK and BK channels are calculated in 
the coupling channel model without fitting parameters. The strong mass shifts down for + and 1 + states have 
been obtained, while 1 + and 2 + states remain almost in situ. The masses of + and 1 + states of B s mesons 
have been predicted. 



After the experimental discovery of the 
£> s (2317) and £>,5(2460) mesons [I], a necessity to 
study the chiral dynamics in heavy-light mesons 
became quite clear. The masses of these states 
proved to be much lower than expected values in 
ordinary quark models while their widths were 
surprisingly small. The problem was studied in 
different approaches: in relativistic quark model 
calculations on the lattice [5J, in QCD 

Sum Rules |6l7j . in chiral models [819] (for re- 
views see also jlOlllj ) . The masses of D s (0 + ) and 
D s (l + ) in closed-channel approximation typi- 
cally exceed by ~ 140 and 90 MeV their experi- 
mental numbers. The main theoretical goal seems 
for us to understand dynamical mechanism re- 
sponsible for such large mass shifts of the + and 
1 + levels and explain why the position of other 
two levels remains practically unchanged. The 
importance of second fact has been underlined by 
S.Godfrey in [3]. 

The mass shifts of the D s (0 + , 1 + ) mesons have 
already been considered in a number of papers 
with the use of unitarized coupled-channel model 
[12] , in nonrelativistic Cornell model [13] , in semi- 
relativistic model with inverse heavy quark mass 
expansion [14] , and in different chiral models [15]- 
[17] . Here we address again this problem with 
the aim to calculate also the mass shifts of the 
D s (l + ) and P s (0 + , 1 + ) states and the widths of 
the 2+ and 1+ states, following the approach de- 
veloped in [IB] , for which strong coupling to the 
S-wave decay channel, containing a pseudoscalar 
(P) Nambu-Goldstone (NG) meson, is crucially 
important. Therefore in this approach princi- 
pal difference exists between vector-vector (VV) 



and VP (or PP) channels. This analysis of two- 
channel system is performed with the use of the 
chiral quark-pion Lagrangian which has been de- 
rived directly from the QCD Lagrangian [18] in 
the frame of the Field Correlator Method (FCM) 
and does not contain fitting parameters, so that 
the shift of the D*(0+) state ~ 140 MeV is only 
determined by the conventional decay constant 
fx. 

From the common point of view, due to spin- 
orbit and tensor interactions the P-wave multi- 
plet of a HL meson is split into four levels with 
J p = 0+, l£, l£, 2+ PI]. Here we use the nota- 
tion H(L) for the higher (lower) 1 + state because 
a priori one cannot say which of them mostly con- 
sists of the light quark j = 1/2 contribution. In 
fact, starting with the Dirac's P-wave levels, one 
has the states with j — 1/2 and j — 3/2. And 
l£ H eigenstates can be parameterized by intro- 
ducing the mixing angle 4>: 



|i+> = cos0|j = |; 



sm^lj = -), 



1 3 

1 1+) = - sin (f>\ j = - ) + cos <f>\j = - ) . 



(1) 



In the heavy-quark (HQ) limit the states with 
j = | and j = h are not mixed, but for finite tuq 
they can be mixed even in closed-channel approx- 
imation [8TT9] . 

Taking the meson emission to the lowest order, 
one obtains the effective quark-pion Lagrangian 
in the form 



AL 



FCM = 



V\ + (x) CT | X 1 75 f^. fa [x ) d 4 X . 



(2) 
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Writing the equation © as AL FC m = — JVifdt, 
one obtains the operator matrix element for the 
transition from the light quark state i (i.e. the 
initial state i of a HL meson) to the continuum 
state / with the emission of a NG meson ((^ a A a ). 
Thus we are now able to write the coupled chan- 
nel equations, connecting any state of a HL me- 
son to a decay channel which contains another 
HL meson plus a NG meson. 

Consider a complete set of the states |/) in the 
decay channel 2 and the set of unperturbed states 
|z) in channel 1. One arrives at the nonlinear 
equation for the shifted mass: 



m\i\ 



m (0) \i] 



E 



i<ot>i- 

Er — m\i] 



(3) 



where to' ' [i] is the initial mass, calculated in 
the single-channel approximation (assumed to be 
known), m[i] - is the final one, Ef is the energy 
of the final state, and the operator V provides the 
transitions between the channels. Note, that in 
our approximation we do not take into account 
the final state interaction in the DK system and 
neglect the D-meson motion. Also, in the w.f. we 
neglect possible (very small) mixing between the 
D{l~i 2 ), D{l~i 2 ) states and between D s (2^ 2 ), 

D 8 (2^, 2 ) states; physical D s (l + ) states can be 
mixed, though. 

In subsequent analysis it is convenient to de- 
fine the masses with respect to nearby threshold: 
fithr = Tor- + Tod. So, we introduce the following 
notations: 

E = to (0) [D s ] 

(4) 

A = Eq + Sm = m[D s \ — mu — mx, (5) 

where A determines the deviation of the D s me- 
son mass from the threshold, and can be complex 
if a decay to DK pair is allowed. In what follows 
we consider unperturbed masses too(J p ) of the 
(Qq) levels as given (our results do not change 
if we slightly vary their position, in this way the 
analysis is actually model- independent). 

For further calculations we should insert the 
explicit meson w.f. to the matrix element in ([3]). 
In our approximation for a HL meson we consider 



a light q (or strange s) quark with current (pole) 
mass to 9 , s moving in the static field of a heavy an- 
tiquark Q, and take its w.f. as a 4-spinor obeying 
the Dirac equation with the linear scalar poten- 
tial and the vector Coulomb potential with frozen 
a s — const: 



U = err, Vc = — 



4 



(6) 



Finally, after long cumbersome calculations 
which are omitted here, we arrive at the the fol- 
lowing equations to determine meson masses and 
widths: 



C s (0 4 



£ [0 + ]-A = j£o(A), 



D s (lt) : 

E [l+] - A = cos 2 <j} ■ T {A) + sin 2 • £2 (A), 

E [l+] - A = sin 2 • ^o(A) + cos 2 <f> ■ ^ 2 (A), 

= sin 2 • f (A) + cos 2 (f> ■ f 2 (A), 
D s (2+ 2 ) : 

E [2+ 2 ] - A = I • jr 2 (A), r[2+ /2 ] = | • f 2 (A), 

where J^o,2 and To j2 are some universal functions; 
definition of those, together with calculation de- 
tails, can be found in [20] . 

In our analysis the 4-component (Dirac) struc- 
ture of the light quark w.f. is crucially important. 
Specifically, the emission of a NG meson is accom- 
panied with the 75 factor which permutes higher 
and lower components of the Dirac bispinors. For 
the j = 1/2, P -wave and the j = 1/2, S -wave 
states it is exactly the case that this "permuted 
overlap" of the w.f. is maximal because the lower 
component of the first state is similar to the 
higher component of the second state and vice- 
versa, while for the analogous overlap between 
j = 3/2, P -wave and the j = 1/2, S -wave states 
the situation is opposite. In the end, it leads to 
the functions J-q, Tq being much larger than Ti, 
r 2 for almost all reasonable values of A. Thus 
the large shift of the 1 + state with a concurrent 
small one for 1 + state reveals a natural explana- 
tion (see below). 
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Now we turn directly to the mass computa- 
tions. We will take into account the following 
pairs of mesons in coupled channels (i refers to 
first (initial) channel, while / refers to second (de- 
cay) one): 



D s (0+) 


D(0-)- 


hif(o-) 


D s (l + ) 




■\-K(0~) 


D s (2+) 


D*(l~) 





(8) 



and analogously for B-meson case, with corre- 
sponding masses and threshold values (in MeV): 

m D + — 1869, m D + + m^- = 2363, 

m D ,+ = 2010, m D ,+ +m K - = 2504, 

m B + = 5279, m B + + m K - — 5772, 

tub* = 5325, tub* + m K - — 5819. 

The light quark eigenfunction is calculated nu- 
merically via Dirac equation with the following 
set of parameters: (the same as in our previous 
papers [21]): 

<7 = 0.18 GeV 2 , a s = 0.39, , , 

10 

m s = 210 MeV, m, - MeV, 

The choice of a and a s is a common one in the 
frame of the FCM approach, and the value of 
the light quark mass really does not influence 
here on any physical results because of its small- 
ness in comparison with the natural mass scale 
\fa. The strange quark mass is taken from [22], 
where it was found from the ratio of experimen- 
tally measured decay constants f(D s )/f(D); the 
same value can be obtained by a renormalization 
group evolution starting from the conventional 
value m s {2 GeV) = 90 ± 15 GeV. 

The ultimate results of our calculations are pre- 
sented in Tables Q3 [3] A priori one cannot say 
whether the \j = \) and \j = §) states are mixed 
or not. If there is no mixing at all, in this case 
the width r(£> a i(2536)) = 0.3 MeV is obtained in 
23J, while the experimental limit is T < 2.3 MeV 
[24] and recently in [25] the width T = 1.0 ± 0.17 
MeV has been measured. Therefore small mixing 
is not excluded and here we take the mixing an- 
gle <j) slightly deviated from (f> = 0° (no mixing 



Table 1 

D s (0 + )-meson mass shift due to the DK decay 
channel and _B s (0 + )-meson mass shift due to the 
BK decay channel (all in MeV) 

state rnW TO ( thc ° r ) m ( c *p) 6m~ 

£> s (0+) 2475 (30) 2330(20) 2317 -145 
B s (0+) 5814(15) 5709 (15) not seen -105" 



case). Then we define those angles cf> which are 
compatible with experimental data for the masses 
and widths of both 1 + states. 

The large value cos 2 <j> for the Vtj(j = 1/2) state 
provides large mass shift (~ 100 MeV) of this 
level and at the same time does not produce the 
mass shift of the lj level, which is almost pure 
j = | state. We would like to stress here that the 
mass shifts weakly differ for D s and B s , or, in 
other words, weakly depend on the heavy quark 
mass. 

Thus we have obtained the shifted masses 
M(B s ,0+) = 5710(15) MeV and M(B s ,l+') = 
5730(15) MeV, which are in agreement with the 
predictions in [TT] and of S.Narison [7] and by 
- 100 MeV lower than in 0,0. The masses of 
the 2+ and 1 + states precisely agree with experi- 
ment. 
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